Density of states and reflectionless tunneling in NS junction with a barrier 
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The effect of the barrier on the proximity effect in normal-superconductor junction is analyzed. A 
general criterion for the barrier, though large, to be effectively transparent, is given. This criterion 
is applied to both the conductance of a disordered NS junction (reflectionless tunneling) and the 
density of states across it, showing that both phenomena stem from the same physical effect. 

PACS numbers: 74.50.-Fr, 74.80.Fp, 73.40.-c 



(N 

"a 
o 
o 



X3 

o 

> 
oo 

0^ 
O 

On 
-I— > 



I 

O 

o 



X 
S3 



At the interface between a superconductor and a nor- 
mal metal Cooper pairs in the superconductor are being 
transformed to electrons and holes in the normal metal 
by the process of Andreev reflection Q. This causes the 
pair amplitude {ipiip]') to be finite in the normal metal 
close to the superconductor (the proximity effect). The 
existence of a barrier at the interface reduces the prob- 
ability of Andreev reflection, and therefore weakens the 
proximity effect. 

The effect of the barrier on the conductance of an NS 
junction (reflectionless tunneling) and on the density of 
states (DOS) across the junction was studied intensively 
(see e.g. |^|-^). It was found that the condition for the 
barrier to substantially reduce the conductance of the 
system is T <C In/d, where F is the transmission proba- 
bility through the barrier, In is the mean free path and d 
is the length of the normal metal between the interfaces 
with the superconductor and the reservoir, van Wees 
et al. used a path integral picture to explain this ef- 
fect. They show that when the disorder in the normal 
metal generates multiple reflections from the interface, 
the constructive quantum interference which results from 
the phase conjugation between electrons and holes en- 
hances the probability of Andreev reflection. We follow 
this treatment and show that the general criterion for the 
barrier not to be effective is that trajectories that reflect 
from the interface more than times before electron- 
hole coherence is lost are dominant. We then show how 
this criterion is applied to the different cases of reflection- 
less tunneling and local DOS with and without magnetic 
field, to give the specific conditions for the barrier to be 
effective in each. We show that in all the different situ- 
ations this criterion reduces to the condition F ^ ln/£., 
where ^ is a characteristic length which is energy, e, and 
magnetic field, H, dependent, and defines the distance 
from the interface, at which electrons and holes are still 
coherent. For the problem of the DOS this gives the same 
condition obtained by Volkov j|] and later by Zhou et al. 
1^. This criterion can be applicable to other systems as 
well, including special cases of ballistic systems as the 
one considered by Hekking and Nazarov |^ . 

In the next section we give a simple and complete treat- 
ment of the effect of the barrier on the DOS in an NS 



junction, using quasiclassical Green function formalism. 
Then we use path integral picture and obtain analytically 
the conductance in the limit of zero temperature, volt- 
age and magnetic field. Finally, we use the path integral 
picture to explain all the above phenomena on the same 
grounds. 

Density of States in NS structure with a barrier. The 
model we consider is an NS slab of width < $s in the 
y direction (^^ = ^JhlXjK is the superconducting coher- 
ence length, Ds is the superconductor diffusion constant, 
A is the superconductor order parameter). The normal 
metal and the superconductor are semi infinite in the x 
direction (perpendicular to the interface). We restrict 
ourselves to the dirty limit, fij/S. > Tn,Ts, where r„ and 
Ts are the normal state elastic mean free times of the 
normal metal and the superconductor. 

Since we are considering the dirty limit we use the 
Usadel equation |9|. We neglect inelastic scattering and 
scattering from paramagnetic impurities and use the sim- 
plest model without self-consistency taking A to be con- 
stant in the superconductor and zero in the normal metal. 
We consider small magnetic field, W Ib {Ib '^s the mag- 
netic length), in which case the derivative of O in the 
y direction can be neglected, and after integrating the 
equation in the y direction we obtain a one dimensional 
equation for B(e, x) [|l0| 

-;Je.^ sine cose + i7sin2e, (1) 

where e(e, x) is defined by the retarded Green functions: 
Ge = cose,F, = sine, and 7 = = W^/M%. 

After solving the equation, the local DOS is obtained 
from N{e, x) = NoReGe{x), Nq is the normal metal DOS 
at the Fermi energy. 

We consider first the case H = 0. Applying the bound- 
ary conditions at -l-oo (normal metal) and —00 (super- 
conductor) one obtains a general solution to Eq. ([T]) px| ] 

4tan"i{tan(^)exp[-a:/C„(e)]}, a; > 



e(e,x) 



e + 4tan-i{tan(^)exp[2:/G(e)]}, x < 0. 



(2) 



Here C„ = y/i/2^n, where ^„ = ^/hDnje is the coher- 
ence length of electron-hole pairs in the normal metal. 
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which determines the energy dependent length scale of 
the proximity effect in the normal metal, and diverges at 
small energies (Dn is the diffusion constant in the normal 
metal). = Cs(l — e^/A^)"^/"* is the decay length for 
states with e < A in the superconductor. The constants 
ipn and "08 are to be determined by the two boundary 
conditions at the NS interface, and Q — tan(A/ — ie) is 
the value of in a bulk superconductor. 

One of the boundary conditions for the Usadel equa- 
tion is a consequence of current conservation, and 
is therefore independent of the existence and strength 
of the barrier at the interface. This boundary condition 
reads: asd/dx{Q{ejO-)) = and/dx{Q{e,0+)), where a„ 
and as are the normal state conductances of the N and S 
metals. The second boundary condition depends on the 
transmission probability of the barrier, and in the limit 
of small transmission reads 11121: 



/„ — e(e,0+) = f sin(e(e,0+) - e(e,0_)). (3) 
ax 

Here F is a measure of the transparency of the barrier, 

and is given by F = | (cos (y9F;^/(l — F^)), where is 

the angular dependent transmission probability of the 

barrier, ip is the angle of incidence at the barrier, and the 

averaging is over the total solid angle. Using the above 

boundary conditions and defining: 77 ~ crn^s/cs^m P — 

,1/2 



for ifjs and ?A„: 



, we obtain two equations 



2/3sin(V'„/2) = F sin(0„ ~ V. ~ 6), (4) 



sin(V's/2) = -77asin('0„/2). (5) 

From Eq. and the definitions of /3 and C„ we see that 
the important parameter of the problem is the ratio be- 
tween F and the parameter In/^n measuring the "dirti- 
ness" of the normal metal 

In the limit of infinite barrier (F/|/3| — 0), ipn — 0, and 
using Eq. (||) ips = 0, and the DOS on each side of the in- 
terface equals its bulk value (zero for the superconductor, 
iVo for the normal metal, see Fig. |l[ thick line). When 
the transmission of the barrier is small compared to the 
"dirtiness" (F/|/3| <C 1), ipn and V's are small, we obtain 
a large jump in the DOS across the barrier, and the de- 
viation from the infinite barrier result is of first order in 
F/|/3| in the superconducting side, and of second order in 
f in the normal metal side (Fig. |l|, dotted line). In the 
limit of |/3|/f = Eq. (|) takes the form of = + 
which is the boundary condition for the case of no barrier 
at the boundary, and since the second boundary condi- 
tion does not depend on the barrier, the results for the 
local DOS in this case are equivalent to those obtained 
for the case of no barrier |jl^,^, even though F <C 1 
(Fig. |l|, thin continuous line). For |/3|/f ^ 1 we have 
tpn — i's — Q ^ 1, and the local DOS is close to the 



no barrier limit (Fig. |T], dashed line). This means that 
in this limit the barrier, though large, has only a small 
effect on the DOS in the system. Since ^„ is energy de- 
pendent and diverges at low energies, there is always an 
energy region where the limit F > Z„/^„ is satisfied and 
the barrier is not effective. In the case F > ln/(,s the 
barrier is not effective for all energies smaller than A. 



FIG. 1. Density of states as a function of distance from the 
NS boundary in units of (negative coordinate values corre- 
spond to the superconductor).e — 0.2 A in all curves. Different 
curves correspond to different ratios of r/|/3|. Thick line - 0, 
dotted - 0.4, dashed - 5, thin line - 00 (no barrier), rj — 1 in 
all curves. 

In the next sections we show that the barrier is not 
effective in certain limits due to the coherence of elec- 
trons and holes in the normal metal. Therefore, it is 
expected that an applied magnetic field will destroy the 
effect. We first consider the effect of finite magnetic 
field at the Fermi energy. In the normal metal the Us- 
adel equation can be solved exactly, and the solution is 
O = 2tan~^[tan(7/'n/2) exp(— ^/7a:)]. In order to solve 
the Usadel equation analytically in the superconductor 
we consider small magnetic fields, H ^ ^q/W^s i-e. 

^ where <I>o is the flux quantum and is defined 
by HW^H = ^0- We then use the boundary condition 
(^, and consider f ^ In/^s- We see that the barrier is 
not effective for F ^ ln/£,H, and for F ^ In/^H the effect 
is destroyed (the barrier becomes effective) Q . 

With both finite magnetic field and finite energy the 
Usadel equation is not solved exactly in the normal metal. 
However, multiplying the equation by dO/dx, integrat- 
ing over X and using the boundary condition at infinity, 
we obtain the first order equation: 

^ - - [^(1 - cose) + 17(1 - cos 26)] i (6) 

Using the boundary condition (^ we obtain an equation 
for the values of O at the interface on both sides: 

— 4ie 1 1 

in[-^{i - cose„(o)) + -7(1 - COs2e„(0))]5 

= fsin(e,(0)-e„(0)). (7) 

To obtain analytical results we assume e <C A and 
H <C ^o/W^s, and then it can be shown that 63(0) w 
7r/2. The condition for the barrier to be noneffective is: 
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f > ln{^-* + G")^''"- This condition reduces to the 
conditions found above for the case of no magnetic field 
and the case of finite magnetic field at zero energy. We 
find that in magnetic field and energy space there is a 
line given by (e/e,)^ + (H/H,)^ = 1, (e, = D^tyil 
He = T^o/lnW) that separates the domain where the 
barrier is not effective (small e and H) from the domain 
where the energy and magnetic field destroy the effect, 
and the barrier is effective. 

Reflectionless tunneling. The phenomenon of reflec- 
tionless tunneling, related to the experimental observa- 
tion of excess low voltage conductance in NS junctions 
(see, e.g. jr^ , p^ ), is treated in many theoretical works 
(e.g. A physically transparent approach to reflec- 

tionless tunneling was given by van Wees et al. . They 
discuss the phenomenon of reflectionless tunneling in a 
semiclassical model, in which the motion of the electrons 
and holes in the normal metal follows a deterministic tra- 
jectory which depends on the initial position and direc- 
tion of the electron, but the barrier between the normal 
metal and the superconductor is treated quantum me- 
chanically. The validity of the semiclassical treatment is 
discussed both in their paper and by Beenakker and van 
Houten |§. 

van Wees et al. consider an N'NS system with a barrier 
at the NS interface. N' is an electron reservoir. Assum- 
ing, at first, a totally reflecting barrier at the interface 
with the superconductor, an electron leaving the reser- 
voir into the disordered normal metal at a given position 
and direction has a given trajectory in the normal metal 
till it will reach the reservoir again. This trajectory can 
not hit the barrier at all or hit the barrier any number of 
times N (see Fig. |). 



superconductor 




reservoir 



FIG. 2. Geometry of the model, an example of a trajectory 
with N=2 

For a disordered normal metal of length d the trans- 
mission probability of electrons across it is T = In/d, and 
therefore, the fraction of trajectories which return to the 
reservoir after hitting the barrier N times is: 



F{N) = 



1 -T 



0, 



(8) 



Considering a barrier with finite transmission probability 



r, at each reflection from the interface either normal or 
Andreev [Q reflections are possible. 

When e — 0, H ~ 0, the Andreev reflected hole retraces 
exactly the path of the incoming electron. Therefore, an 
incoming electron can return to the reservoir either as an 
electron at the end of the trajectory, or as a hole at the 
beginning of the trajectory (Fig. ^). 

We define ree{N) and rhh{N)) as the amplitudes of 
electron and hole reflections in an N trajectory, re- 
spectively. The average contribution to the charge cur- 
rent of trajectories involving N reflections at the inter- 
face is I{N) — 2\rhe{N)\'^ and therefore the conduc- 
tance of the system is given by: G{V ^ Q,H — Q) — 
X]w=o -^(-^)^(^)' where Gs is the Sharvin conduc- 
tance of the interface, and for a system with a finite 
cross-section Gs = 2e^n/h, n is the number of trans- 
verse channels. Considering first the case where e = 0, 
-ff 0, we obtain a recurrence formula: rhe{N) ~ 
i{\rhe{N-l)\ + \rhe\)/{l + \rhe\\rhe{N where rhe 
is the amplitude for the Andreev reflection of an electron 
at the interface, given by The = r^h = iF/(2 — F). The 
above formula is obtained by summing the amplitudes of 
all the different paths within the trajectory that result 
in a hole being reflected back to the reservoir (a path is 
a given sequence of normal and Andreev reflections from 
the different points the trajectory intersects the interface. 
There are infinitely many paths for a given trajectory). 

We then obtain: rhe{N) = i tanh(Af tanh"\|r/,.e(l)|)) 
and for a barrier with small transmission probability 

rhe{N) « ^ tanh(A^|r,,e(l)|) « i tanh(7VF/2). (9) 

Inserting this result in the conductance formula we ob- 
tain: 



G{V ^ 0,i/ = 0) 



4^(^ + i)-i (F»r). ^ 



(r « T), 



This result differs from the result obtained by Beenakker 
et al. 1^ by a factor of 2 in the F 3> T limit since they con- 
sider the case of short range disorder while we consider 
the case where the electrons and holes follow a classical 
trajectory, which prevails when the scattering potential 
varies slowly on the scale of a wavelength [p^ . 

The strong dependence of the amplitude of Andreev 
reflection on the number of times the trajectory hits the 
interface, as manifested in Eq. (^), is due to the coher- 
ence of the incoming electron and reflected hole. Due 
to this coherence, the probability of Andreev reflection, 
\rhe{N)\ approaches 1 for N 3> F^^, and the probability 
of normal reflection, |ree(A^)| approaches zero. The phe- 
nomenon of reflectionless tunneling is a result of the fact 
that as the disorder is enhanced, the fraction of large N 
trajectories grows. 

Reflectionless tunneling and density of states. Both 
reflectionless tunneling and the influence of the barrier 
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on the DOS in an NS structure are related to the co- 
herent transport of pairs from the superconductor across 
the barrier to the normal metal, which is given by the 
probability of an electron approaching the boundary to 
be Andreev reflected from it. The probability of Andreev 
reflection approaches 1 for trajectories with N ^ F^^, 
so the barrier is not effective if large N trajectories occur 
with high probability. Therefore, the general criterion 
for the barrier not to be effective is that trajectories that 
hit the barrier more than F"^ times before losing the 
coherence between holes and electrons are dominant in 
the system. Using random walk theory it can be shown 
that the length of a trajectory between N consecutive 
times it hits the barrier is of the order of Zat « 7V^Z„ 
1^^ , and therefore only when coherent trajectories with 
lengths larger than Lr = In/^'^ occur with high proba- 
bility will the barrier not be effective. This criterion can 
be applied to various cases, where different mechanisms 
limit the length of coherent trajectories. 

In the case of reflectionless tunneling the length of 
the normal metal, d, is what limits the trajectories to 
lengths of order <P /In and therefore the barrier is not 
effective when F ^ T. In the case of the DOS in 
NS structure with semi-infinite normal metal, there is 
no limit to the length of the trajectories, and there- 
fore, when e = and H = 0, the barrier is not effec- 
tive. However, at finite e and H the electron-hole coher- 
ence is limited. The equation for r/ie(2) is generalized 
g: rhe{2) = [rfte(l + exp(iA0)]/[l + |r;,e||reh|exp(iA(^)]. 
Here Acj) = 2eL/?ivp+4:Tr{HA/^o), the first part is due to 
the phase accumulated by an electron and a hole travers- 
ing a part of a trajectory of length L between the two 
intersection points with the barrier, and the second part 
is due to the magnetic flux, HA, through the Andreev 
loop consisting of the part of the trajectory of length 
L, and closed by the superconductor. For electrons at 
e = 0,i7 = 0, = 0, rhe{2) = 2rfte/(l + Ir/^eP) and 
the probability I r/ie (2) p has a maximum. In the same 
way \rhe{N)\'^ (g) has a maximum when Acj) ~ 0. When 
electron-hole coherence breaks (Acj) sa 2tt) there is no en- 
hancement in the Andreev reflection probability due to 
constructive interference, and the barrier becomes effec- 
tive. At flnite energies coherence will be maintained up 
to lengths of order Lc = Tivp/e, and therefore the DOS 
only at e <C F^?i/T„ will not be affected by the barrier 
(Lr ^ Lc). This is the same condition (F ^ ln/£,n) 
which was obtained using the Usadel equation. 

At e = 0, 7^ 0, it can be shown that in a slab of 
width W a trajectory of length 1%/Wln encloses one flux 
quantum and therefore the condition for the barrier not 
to be effective is F 3> Uil^Hi again as was obtained using 
the Usadel equation. 

Similar considerations lead to the result that reflec- 



tionless tunneling is destroyed at energies of the order of 
the Thouless energy, hDn/(P, and magnetic flelds of the 
order oi H = ^o/Wd. 

The criterion of hitting the barrier more than F^^ 
times before losing the electron-hole coherence can also 
be applied to ballistic systems where the geometry al- 
lows returns to the interface, as an SININ double barrier 
system |18|, or the system considered by Hekking and 
Nazarov |p|. 
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